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The review of current experimental situation in the measurements of the first 
moment T p ^ n of spin dependent nucleon structure functions gi- Pjn (x, Q 2 ) is presented. 
The results of the calculations of twist-4 corrections to r Pi „ are discussed and their 
accuracy is estimated. The part of the proton spin X carried by u, d, s quarks 
is calculated in the framework of the QCD sum rules in the external fields. The 
operators up to dimension 9 are accounted. An important contribution comes from 
the operator of dimension 3, which in the limit of massless u, d, s quarks is equal 
'-^ ' to the derivative of QCD topological susceptibility The comparison with 

the experimental data on £ gives x'(0) = (2.3 ± 0.6) x 1CT 3 GeV 2 . The limits on 
K> . £ and x'(0) are found from selfconsistency of the sum rule, £ > 0.05, x'(0) ^ 

1.6 x 10~ 3 GeV 2 . The values of g A = 1.37 ± 0.10 and g\ = 0.65 ± 0.15 are also 
determined from the corresponding sum rules. 
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/ dedicate this lecture to the memory of my friend Volodya Gribov, whom I knew for 
about half a centure. Now it becomes even more clear how great was his influence on 
physics: his brilliant ideas, his uncompromising approach to science, his teaching ability. 
My loss is even more painful: every meeting with Volodya was like a holyday to my soul. 

1. Introduction. Recent experimental data. 

In the last years, the problem of nucleon spin content and particularly the question 
which part of the nucleon spin is carried by quarks, attracts a strong interest. The valu- 
able information comes from the measurements of the spin-dependent nucleon structure 
functions g\{x, Q 2 ) in deep inelastic e(fi)N scattering (for the recent data see [1,2,3], for a 
reviews [4,5]). The parts of the nucleon spin carried by u, d and s-quarks are determined 
from the measurements of the first moment of gi(x, Q 2 ) 

i 

r p , n (<5 2 ) = J dxg 1;Pin (x,Q 2 ) (1) 



At high Q 2 with the account of twist-4 contributions T P:fl (Q 2 ) have the form 

r P AQ 2 ) = KUQ 2 ) + K w n(Q 2 ) ( 2 ) 



^;„(Q 2 ) = ^{[i 



a - 3.58a 2 - 20.2a 3 - ca 4 ][±g A + \g\] 



4 1 I N* 

+ -[l--a- 0.55a 2 - 4.45a 3 ]E j - -^-a s (Q 2 )Ag(Q 2 ) (3) 

In eq.(3) a = a s (Q 2 )/n, g A is the /3-decay axial coupling constant, g A = 1.260 ± 0.002 [6] 

g A = Au - Ad g\ = Au + Ad - 2 As E = Au + Ad + As. (5) 

A-u, Ad, As, Ag are parts of the nucleon spin projections carried by u, d, s quarks and 
gluons: 



Aq 



q + (x) - q-{x) 



dx (6) 



where q+(x), q~(x) are quark distributions with spin projection parallel (antiparallel) to 
nucleon spin and a similar definition takes place for Ag. The coefficients of perturbative 
series were calculated in [7-10], the numerical values in (3) correspond to the number of 
flavours Nf = 3, the coefficient c was estimated in [11], c « 130. In the MS renormal- 
ization scheme chosen in [7-10] g A , g A and E are Q 2 -independent. In the assumption of 
the exact SU (3) flavour symmetry of the octet axial current matrix elements over baryon 
octet states g A = 3F — D = 0.59 ± 0.02 [12]. On the basis of operator product expansion 
(OPE) the quantities g A ,g A and E are related to the proton matrix element of isovector, 
octet and singlet axial currents correspondingly: 
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2ms^g A , g\, £) = (P, s | \p,s), (7) 

where is the proton spin 4-vector, m is the proton mass. 

Strictly speaking, in (3) the separation of terms proportional to E and Ag is arbitrary, 
since OPE has only one singlet in flavour twist-2 operator for the first moment of the polar- 
ized structure function - the operator of singlet axial current jj^ix) = Qi( x )l^l5<l, Q = 

i 

u, d, s. The separation of terms proportional to E and Ag is outside the framework of 
OPE and depends on the infrared cut-off. The expression used in (3) is based on the 
physical assumption that the virtualities p 2 of gluons in the nucleon are much larger than 
light quark mass squares, \p 2 \ ^> m 2 [13] and that the infrared cut-off is chosen in a way 
providing the standard form of axial anomaly [14]. 

Since the separation from E of the term, proportional to Ag, results in redefinition 
of E, sometimes in the analysis of the data it is separated, sometimes it is not. In what 
follows in the main part of the Lecture I will not separate Ag contribution from E, only 
sometimes mentioning how large it could be. 

Twist-4 corrections to T pn were calculated by Balitsky, Braun and Koleshichenko 
(BBK) [15] using the QCD sum rule method. 

BBK calculations were critically analyzed in [16], where it was shown that there are 
many possible uncertainties in these calculations: 1) the main contribution to QCD sum 
rules comes from the last accounted term in OPE - the operator of dimension 8; 2) there 
is a large background term and a much stronger influence of the continuum threshold 
comparing with usual QCD sum rules; 3) in the singlet case, when determining the induced 
by external field vacuum condensates, the corresponding sum rule was saturated by r\- 
meson, what is wrong. The next order term - the contribution of the dimension 10 
operator to the BBK sum rules was estimated by Oganesian [17]. The account of the 
dimension-10 contribution to the BBK sum rules and estimation of other uncertainties 
results in (see [16]): 

b p _ n = -0.006 ± 0.012 GeV 2 (8) 

b p+n = -0.035(±100%) GeV 2 (9) 

As is seen from (8), in the nonsinglet case the twist-4 correction is small (< 2% at 
Q 2 > 5GeV 2 ) even with the account of the error. In the singlet case the situation is much 
worse: the estimate (9) may be considered only as correct by the order of magnitude. 

One may expect that at low Q 2 < 3 GeV 2 the nonperturbative (higher twist) cor- 
rections to T p ^ n (Q 2 ) are much larger in absolute values, than given by (8), (9). This 
statement follows from the requirement, that at Q 2 = T p ^ n (Q 2 ) satisfies the Gerasimov- 
Drell-Hearn (GDH) sum rule and a smooth connection of T P:n (Q 2 ) at intermediate Q 2 
and those at Q 2 = should exist. (In accord with the GDH sum rule r p , n (0) = and 
TpnW = -K pn/^ m2 » where K PtJl are proton and neutron anomalous magnetic moments 
- see [16].) In [16] the model was suggested, which realizes such smooth connection. As 
was demonstrated in [16] the model is in a good agreement with the recent experimental 
data. An interesting feature of the model, supported by the data, is that the sign of 
nonperturbetive correction coincides with the sign of twist-4 terms (7), (8) in the case of 
proton, but it is opposite for neutron. 
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I turn now to comparison of the theory with the recent experimental data. In Table 
1 the recent data obtained by SMC [1], E154(SLAC) [2] and HERMES [3] groups are 
presented. 

Table 1 





r P 


r 

L n 


r — r 


a s (5GeV 2 ) 


SMC 


0.132 ±0.017 


-0.048 ± 0.022 


0.181 ±0.035 


0.270^ 


combined 


0.142 ±0.011 


-0.061 ±0.016 


0.202 ±0.022 


0.116^- 


E 154(SLAC) 


0.112 ±0.014 


-0.056 ±0.008 


0.168 ±0.012 


0.339 iffi 


HERMES 




-0.037 ±0.015 






EJ/Bj sum rules 


0.168 ±0.005 


-0.013 ±0.005 


0.181 ±0.002 


0.276 



In the second line of Table 1 the results of the performed by SMC [1] combined analysis 
of SMC [1], SLAC-E80/130 [18], EMC [19] and SLAC-E143 [20] data are given. The data 
presented in the first three lines of Table 1 refer to Q 2 = 5 GeV 2 , HERMES data refer 
to Q 2 = 2.5 GeV 2 . In all measurements each range of x corresponds to each own mean 
Q 2 . Therefore, in order to obtain gi(x, Q 2 ) at fixed Q 2 the authors of ref.'s [1,2] used the 
following procedure. At some reference scale Ql (Ql = lGeV 2 in [1] and Ql = 0.3AGeV 2 
in [2]) quark and gluon distribution were parametrized as functions of x. (The number 
of the parameters was 12 in [1] and 8 in [2]). Then NLO evolution equations were solved 
and the values of the parameters were determined from the best fit at all data points. 
The numerical values presented in Table 1 correspond to MS regularization scheme, 
statistical, systematical, as well as theoretical errors arising from uncertainty of a s in the 
evolution equations, are added in quadratures. The HERMES value of T n , measured at 
Q 2 = 2.5 GeV 2 can be recalculated to Q 2 = 5 GeV 2 using the model [16], matching GDH 
sum rule at Q 2 = and asymptotic behavior of T n (Q 2 ). The result is: T n (Q 2 = 5 GeV 2 ) = 
-0.045 ± 0.015 (HERMES). In the last line of Table 1 the Ellis- Jaffe (EJ) and Bjorken 
(Bj) sum rules prediction for T p , T n and T — T n , correspondingly are given. The EJ sum 
rule prediction was calculated according to (3), where As = , i.e., £ = g\ — 0.59 was 
put and the last-gluonic term in (3) was omitted. The twist-4 contribution was accounted 
in the Bj sum rule and included into the error in the EJ sum rule. The a s value in the 
EJ and Bj sum rules calculation was chosen as a s (5 GeV 2 ) = 0.276, corresponding to 
a s (M z ) = 0.117 and = 3Q0MeV (in two loops). As is clear from Table 1, the 

data, especially for T n contradict the EJ sum rule. In the last column, the values of a s 
determined from the Bj sum rule are given with the account of twist-4 corrections. 

The experimental data on T p presented in Table 1 are not in a good agreement. 
Particularly, the value of T p given by E154 Collaboration seems to be low: it does not 
agree with the old data presented by SMC [21] (T p = 0.136 ± 0.015) and E143 [20] 
(r p = 0.127±0.011). Even more strong discrepancy is seen in the values of a S) determined 
from the Bj sum rules. The value which follows from the combined analysis is unacceptably 
low: the central point corresponds to = 15MeV\ On the other side, the value, 

determined from the E154 data seems to be high, the corresponding a s (M z ) = 0.126 ± 
0.009. Therefore, I come to a conclusion that at the present level of experimental accuracy 
a s cannot be reliably determined from the Bj sum rule in polarized scattering. 

Table 2 shows the values of £ - the total nucleon spin projection carried by u, d and 
s-quarks found from T p and T p presented in Table 1 using eq.(3). (It was put qa = 
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1.260, g\ = 0.59, the term, proportional to Ag is included into £.). 



Table 2: The values of E 





From T p 


From r n 




At a s (5GeV 2 ) = 
= 0.276 


At a s (5GeV 2 ) 
given in Table 1 


At a s (5GeV 2 ) = 
= 0.276 


At a s (5GeV 2 ) 
given in Table 1 


SMC 


0.296 


0.294 


0.294 


0.296 


Comb. 


0.390 


0.290 


0.175 


0.255 


E154 


0.110(0.17; 0.29) 


0.17(0.24; 0.34) 


0.22(0.28; 0.17) 


0.17 (0.24; 0.13) 


HERMES 






0.38(0.26) 


at a s (2.5GeV 2 ) = 
= 0.337 



In their fitting procedure [2] E154 Collaboration used the values g\ = 0.30 and gA = 
1.09. The values of E obtained from T p and T n given by E154 at g\ = 0.30,^ = 1.26 
and g\ = 0.30, g& = 1.09 are presented in parenthesis. The value g\ = 0.30 corresponds 
to a strong violation of SU(3) flavour symmetry and is unplausible; gA = 1-09 means a 
bad violation of isospin and is unacceptable. As seen from Table 1, E is seriously affected 
by these assumptions. The values of E found from T p and T n using SMC and combined 
analysis data agree with each other onlyif one takes for a s (5GeV 2 ) the values given in 
Table 1 (a s = 0.116 for combined data), what is unacceptable. 

The twist-4 corrections were accounted in the calculations of E in Table using eq.'s 
(8), (9). At Q 2 = 5 GeV 2 they result in increasing of E by 0.04 if determined from T n , at 
Q 2 = 2.5 GeV 2 (HERMES data) the twist-4 correction increase E by 0.06. In the last line 
in parenthesis is given the value of E, when higher twist corrections were found basing 
on the model matching GDH sum rule and asymptotic behavior of r n [16]. The chosen 
value of a s (2.5 GeV 2 )=0.337 corresponds to the same Aq^ d (2 loops) = 360 MeV, as 
a s (5 GeV 2 ) = 0.276. 

To conclude, one may say, that the most probable value of E is E ss 0.3 ± 0.1 . 
The contribution of gluons may be estimated as Ag(lGeV 2 ) 0.3 (see [16]). Then 
Ag(5GeV 2 ) ~ 0.6 and the account of gluonic term in eq.(3) results in increasing of E by 
0.06. At E = 0.3 we have Am = 0.83, Ad = -0.43, As = -0.1. 

2. The QCD sum rules calculation of E. 

The quantity E, which has the meaning of proton spin projection, carried by u, d, s 
quarks is of a special interest. 

An attempt to calculate E using QCD sum rules in external fields was done in ref.[22]. 
Let us shortly recall the idea. The polarization operator 

n(p) = i J d 4 xe tpx (0\T{ri(x),f](0)}\0) (10) 

was considered, where 

T]{x) = e^fu^C^u^y^ix) (11) 

is the current with proton quantum numbers [23], [24] u a ,d b are quark fields, a,b,c are 
colour indeces. It is assumed that the term 

AL = iA (12) 
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where is a constant singlet axial field, is added to QCD Lagrangian. In the weak axial 
field approximation H(p) has the form 

n( P ) = rt \p) + ii£\p)A,. (is) 

nW(p) is calculated in QCD by OPE at p 2 < 0, \p 2 \ 3> R~ 2 , where R c is the confine- 
ment radius. On the other hand, using dispersion relation, U^(p) is represented by the 
contribution of the physical states, the lowest of which is the proton state. The contri- 
bution of excited states is approximated as a continuum and suppressed by the Borel 
transformation. The desired answer is obtained by equalling these two representations. 
This procedure can be applied to any Lorenz structure of H^\p) , but as was argued in 
[25,26], the best accuracy can be obtained by considering the chirality conserving structure 

An essential ingredient of the method is the appearance of induced by the external 
field vacuum expectation values (v.e.v). The most important of them in the problem at 
hand is 

(0\j%\0) A = 3f 2 A, (14) 

of dimension 3. The constant /q is related to QCD topological susceptibility. Using (12), 
we can write 

(0b'»A = Um q ^ i J d 4 xe^(0\T{f u5 (x),j%(0)}\0)A u = 

= lim q ^P^(q)A u (15) 

The general structure of P^ u (q) is 

P,u(q) = -P L (q 2 )S^ + P T {q 2 ){-^ v q 2 + q,q u ) (16) 

Because of anomaly there are no massless states in the spectrum of the singlet polarization 
operator P^ v even for massless quarks. Pt,l(<? 2 ) also have no kinematical singularities at 
q 2 = . Therefore, the nonvanishing value P^ u {fS) comes entirely from Pi(q 2 ). Multiplying 
P^(q) by QtiQv, in the limit of massless u, d, s quarks we get 

q^PM = ~PL(q 2 )q 2 = N 2 {a s /^fi J dV^x 

x (o\tg; u (x)g; u (x),gz(o)gZ(oM, (17) 

where G™ v is the gluonic field strength, G^ = (l/2)e tJLV \ a G\ (T .( r Ih.e anomaly condition was 
used, Nf = 3.). Going to the limit q 2 — > 0, we have 

f 2 = -(1/3)P L (0) = ^Y(O), (18) 
where x(<? 2 ) is the topological susceptibility 

X(q 2 ) = * J d 4 xe^{0\TQ 5 (x) : Q 5 (0)\0) (19) 
and Q^(x) is the topological charge density 
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Q B (x) = (a s /8n) G^(x)G^(x), (20) 

As is well known [27], x(0) = if there is at least one massless quark. The attempt to 
find hself by QCD sum rules failed: it was found [22] that OPE does not converge 
in the domain of characteristic scales for this problem. However, it was possible to derive 
the sum rule, expressing E in terms of /q (14) or x'(0). The OPE up to dimension d — 7 
was performed in ref.[22]. Among the induced by the external field v.e.v.'s besides (14), 
the v.e.v. of the dimension 5 operator 

g(0\J2qia(l/2)\ n G n ^q\0) A ^3h A p , q = u,d,s (21) 
i 

was accounted and the constant h was estimated using a special sum rule, 
h m 3 x 10~ 4 GeV A . There were also accounted the gluonic condensate d = 4 and the 
square of quark condensate d — 6 (both times the external field operator, d — 1). 
However, the accuracy of the calculation was not good enough for reliable calculation of 
£ in terms of /q : the necessary requirement of the method - the weak dependence of the 
result on the Borel parameter was not well satisfied. 

In [28] the accuracy of the calculation was improved by going to higher order terms in 
OPE up to dimension 9 operators. Under the factorization assumption - the saturation of 
the product of four-quark operators by the contribution of an intermediate vacuum state 
- the dimension 8 v.e.v.'s were accounted (times A^): 

- g(0\gaafi(l/2)X n G^q ■ qq\0) = m 2 (0\qq\0) 2 , (22) 

where m\ = 0.8 ± 0.2 GeV 2 was determined in [28]. In the framework of the same 
factorization hypothesis the induced by the external field v.e.v. of dimension 9 

*s(0\j ( 3\0)A(0\qq\0) 2 (23) 

is also accounted. In the calculation the following expression for the quark Green function 
in the constant external axial field was used [26]: 

(0\T{q a a (x), gJ(0)}|0) A = t5 ab x a p/27r 2 x 4 + 
+ (l/27r 2 )5 ab (Ax)( l5 xU/x 4 - (l/12)5 ab 5 aP (0\qq\0) + 



+ (l/72)i5 ab (0\qq\0)(xA l5 - Ax l5 ) af3 + 



+ (1/12) f 2 5 ab (A l5 ) a(S + (1/216)^0 



ab i 



(5/2)rr 2 i 75 - {Ax)x^ 



(24) 



a/3 



The terms of the third power in x-expansion of quark propagator proportional to A^ 
are omitted in (24), because they do not contribute to the tensor structure of n M of 
interest. Quarks are considered to be in the constant external gluonic field and quark 
and gluon QCD equations of motion are exploited (the related formulae are given in [29]). 
There is also an another source of v.e.v. ho to appear besides the x-expansion of quark 
propagator given in eq.(24): the quarks in the condensate absorb the soft gluonic field 
emitted by other quark. A similar situation takes place also in the calculation of the 
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v.e.v. (23) contribution. The accounted diagrams with dimension 9 operators have no 
loop integrations. There are others v.e.v. of dimensions d < 9 particularly containing 
gluonic fields. All of them, however, correspond to at least one loop integration and are 
suppressed by the numerical factor (2it)~ 2 . For this reason they are disregarded. 
The sum rule for E is given by 

E + C M 2 = -1 + -le m2 / M2 {a 2 L A ' 9 + 

W^grnu^gj^-if-i^ (25) 

Here M 2 is the Borel parameter, \ N is defined as X 2 N = 32ir 4 \ 2 N = 2.1 GeV 6 , (0\r)\p) = 
\n v p , where v p is proton spinor, W 2 is the continuum threshold, W 2 = 2.5 GeV 2 , 

a = -(27r) 2 (0|gg|0) = 0.55 GeV 3 (26) 

E (x) = 1 - e~ x , E^x) = 1 - (1 + x)e~ x 

L = ln(M/A)/ln(fx/A), A = Aq C d = 200 MeV and the normalization point /i was 
chosen /i = 1 GeV. 

When deriving (25) the sum rule for the nucleon mass was exploited what results in 
appearance of the first term, -1, in the right hand side (rhs) of (25). This term absorbs 
the contributions of the bare loop, gluonic condensate as well as a s corrections to them 
and essential part of terms, proportional to a 2 and m^a 2 . It must be stressed, that with 
the account of dimension 9 operators the OPE series in the calculation of E is going 
up to the same order as OPE in the calculation of nucleon mass, where in the chirality 
conserving sum rule the operators up to dimension 8 were accounted (see Appendix, one 
additional dimension in the sum rule for E comes from the dimension of external axial 
field A^). Therefore, both sum rules are on the same footing and the procedure of using 
chirality conserving nucleon sum rule (A.l) in (25) is legitimate. Otherwise, and this was 
the drawback of calculations in [25], [26], the approach is not completely selfconsistent. 
The values of the parameters, a, \ 2 N , W 2 taken above were chosen by the best fit of the 
sum rules for the nucleon mass (see [30] and Appendix) performed at A = 200 MeV. 
It can be shown, using the value of the ratio 2m s /(m u + m^) = 24.4 ± 1.5 [31] that 
a(l GeV) = 0.55 GeV 3 corresponds to m s {\ GeV) = 153 MeV. a s corrections are 
accounted in the leading order (LO) what results in appearance of anomalous dimensions. 
Therefore A has the meaning of effective A in LO. Its numerical value does not contradicts 
two loops value of A, used in Sec.l. (Formally, A^ 3 ^(2 loops) = 360 MeV would results to 
Kf) f {LO) = 250 MeV.) 

The unknown constant Cq in the left-hand side (lhs) of (25) corresponds to the contri- 
bution of inelastic transitions p —> N* interaction withA^ — > p (and in inverse order). 
It cannot be determined theoretically and may be found from M 2 dependence of the rhs 
of (25) (for details see [30,32]). The necessary condition of the validity of the sum rule is 
|E| > \C M 2 \exp[(-W 2 + m 2 )/M 2 ] at characteristic values of M 2 [32]. The contribution 
of the last term in the rhs of (25) is negligible. The sum rule (25) as well as the sum rule 
for the nucleon mass is reliable in the interval of the Borel parameter M 2 where the last 
term of OPE is small, less than 10 — 15% of the total and the contribution of continuum 
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does not exceed 40 - 50% . This fixes the interval 0.85 < M 2 < 1.4 GeV 2 . The in- 
dependence of the rhs of (25) at /q = 3 x 10~ 2 GeV 2 is plotted in Fig.l. The complicated 
expression in rhs of (25) is indeed an almost linear function of M 2 in the given interval! 
This fact strongly supports the reliability of the approach. The best values of £ = E^* 
and C = C^ 1 are found from the x 2 fitting procedure 

X 2 = - E W lt + CtM 2 - R(M 2 )] 2 = min, (27) 
n i=i 

where R(M 2 ) is the rhs of (25). 

The values of E as a function of /q are plotted in Fig. 2 together with v5P- I n the used 
above approach the gluonic contribution cannot be separated and is included in E. As 
discussed in Sec. 1 the experimental value of E can be estimated as E = 0.3 ± 0.1. Then 
from Fig.2 we have / 2 = (2.8 ± 0.7) x 10~ 2 GeV 2 and x'(°) = ( 2 - 3 ± °- 6 ) x 10_3 Gey2 
. The error in /q and x\ besides the experimentall error, includes the uncertainty in the 
sum rule estimated as equal to the contribution of the last term in OPE (two last terms 
in Eq.25) and a possible role of NLO a s corrections. At f% < 0.02 GeV 2 x 2 is much worse 
and the fit becomes unstable. This allows us to claim (with some care, however,) that 
x'(0) > 1.6 x 10~ 3 GeV 2 and E > 0.05 from the requirement of selfconsistency of the sum 
rule. The x 2 curve also favours an upper limit for E < 0.6. At /q = 2.8 x 10~ 2 GeV 2 the 
value of the constant C found from the fit is C = 0.19 GeV~ 2 . Therefore, the mentioned 
above necessary condition of the sum rule validity is well satisfied. 

Let us discuss the role of various terms of OPE in the sum rules (25) To analyze it we 
have considered sum rules (25) for 4 different cases, i.e. when we take into consideration: 
a) only contribution of the operators up to d=3 (the term -1 and the term, proportional 
to /q in (25)); b) contribution of the operators up to d=5 (the term ~ h is added); c) 
contribution of the operators up to d=7 (three first terms in (25)), d) our result (25), i.e. 
all operators up to d=9. For this analysis the value of /q = 0.03 GeV 2 was chosen, but the 
conclusion appears to be the same for all more or less reasonable choice of f$. Results of 
the fit of the sum rules are shown in Table 3 for all four cases. The fit is done in the region 
of Borel masses 0.9 < M 2 < 1.3 GeV 2 . In the first column the values of E are shown , in 
the second - values of the parameter C, and in the third - the ratio 7 = |vA^/E|, which 
is the real parameter, describing reliability of the fit. From the table one can see, that 
reliability of the fit monotonously improves with increasing of the number of accounted 
terms of OPE and is quite satisfactory in the case d 



Table 3 



case 


E 


C(GeV-' 2 ) 


7 


a) 


-0.019 


0.31 


10- 1 


b) 


0.031 


0.3 


5.10" 2 


c) 


0.54 


0.094 


9.10" 3 


d) 


0.36 


0.21 


1.3 • 10- 3 



Recently, the first attempt to calculate x'(0) 011 the lattice was performed [33]. The 
result is x'(0) = (0-4 ± 0.2) x 10~ 3 GeV 2 , much below our value. However, as mentioned 
by the authors, the calculation has some drawbacks and the result is preliminary. 
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In the papers by Narison, Shore and Veneziano (NSV) [34], [35], an attempt to find 
the links between E and x'(0) was done. NSV found that £ is proportional to \Jx'(fy 
and calculated °y QCD sum rules. From my point of view, the approach of ref.'s 

[33], [34] is not justifiable. Instead of use of firmly based and self consistent OPE, as was 
done above, in [34], [35] the matrix element < p \ Q 5 \ p > was saturated by contribution of 
two operators Q$ and singlet pseudoscalar operator T>qj 5 q - and the result was obtained 
by orthogonalization of the corresponding matrix. I have doubts that such procedure can 
be grounded. The calculation of x'(0) by QCD sum rules is not correct, because, as was 
shown in [22] by considering in the same problem with account of higher order terms of 
OPE, than it was done in [34], [35], the OPE breaks down at the scales, characteristic for 
this problem. I do not believe, that the value x'(0) = (0.5 ± 0.2) • 10~ 3 GeV 2 found in 
[34] is reliable. 

3. Calculation of proton axial coupling constant g\ and qa- 

From the same sum rule (25) it is possible to find g\ - the proton coupling constant 
with the octet axial current, which enters the QCD formula for T PtTl . There are two 
differences in comparison with (25): 

1. Instead of f{j it appears the square /| of the pseudoscalar meson coupling constant 
with the octet axial current. In the limit of strict SU(3) flavour symmetry it is equal 
to f 2 , f n = 133 MeV. However, it is known, that SU(3) symmetry is violated and the 
kaon decay constant, fx ~ 1.22/^ [6]. In the linear in s-quark mass m s approximation 
f v = 1.28/V. We put for /f the value /f = 2.6 x 10~ 2 GeV 2 , intermediate between f 2 and 
f 2 . 

2. h should be substituted by raff 2 . The constant m\ is determined by the sum rules 
suggested in [36]. A new fit corresponding to the values of the parameters used above, 
was performed and it was found; m\ = 0.16 GeV 2 . 

The M 2 -dependence of g\ + C 8 M 2 is presented in Fig.l and the best fit according to 
the fitting procedure (27) at 1.0 < M 2 < 1.3 GeV 2 gives 

g\ = 0.65 ± 0.15, C 8 = 0.10 GeV 2 yx 2 " = 1-2 x 10 -3 (28) 

(The error includes the uncertainties in the sum rule as well as in the value of /|). The 
obtained value of g\ within the errors coincides with g\ = 0.59 ± 0.02 [12] found from 
the data on baryon octet /9-decays under assumption of strict SU(3) flavour symmetry 
and contradicts the hypothesis of bad violation of SU(3) symmetry in baryon axial octet 
coupling constants [37]. 

A similar sum rule with the account of dimension 9 operators can be derived also for 
gA - the nucleon axial /3-decay coupling constant. It is an extension of the sum rule found 
in [25] and has the form 



g A + C A M 2 = 1 + 



m 2 /M 2 



(29) 



The main term in OPE of dimension 3 proportional to f 2 occasionally was cancelled. For 
this reason the higher order terms of OPE may be more important in the sum rule for 
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g A than in the previous ones. The M 2 dependence of g A — 1 + CaM 2 is plotted in Fig.l, 
lower curve; the curve is almost the straight line, as it should be. The best fit gives 

g A = 1.37 ± 0.10, C A = -0.088 GeV~\ ^ = 1.0 x 10~ 3 (30) 

in comparison with the world average g A = 1.260 ± 0.002 [6]. The inclusion of dimension 
9 operator contribution essentially improves the result: without it g A would be about 1.5 
and x 2 would be much worse. 

The work was supported in part by CRDF Grant RP2-132, INTAS Grant 93-0283, 
RFFR Grant 97-02-16131 and Swiss Grant 7SUPJ048716. 

Appendix 

The fit of the sum rules for nucleon mass. 

Since in comparison with previous fit [30] of the sum rules for nucleon mass the value 
of QCD parameter was changed now, the new fit was performed. (In the previous cal- 
culations it was used A = 100 MeV, now we take A = 200 MeV.). The sum rules for 
chirality conserving and chirality violating parts of the polarization operator II^Hp) (6) 
defined by (3) are correspondingly 

m ^(S) l ~ 4/9 + t a2L4/9 + i 6m2 MS) l ~ 4/9 - 

-\ a2 W 2=l2 ^ m2/M2 (A - 1} 

9nM*F ( W2 \ + 272 as{M2) ° 3 1 nh - m\ 2 P ~^l^ (A 9^ 

where 

b = (2tt) 2 (0 I —G 2 | 0) = 0.50 GeV\ 

7T 

2 

E 2 (x) = 1 - (l + x + ^-)e~ x 

and the other notations are the same as in (25), (26). Parameters a and W 2 were treated 
as fitting parameters and it was required that in the fitting interval 0.8 < M 2 < 1.3 GeV 2 
the quantities X 2 N found from both sum rules (A.l) and (A. 2) must be close to one another 
and close to a constant, independent of M 2 . The values of A^, determined from (A.l) 
and (A. 2) as functions of a (at normalization point ji = 1 GeV and continuum threshold 
W 2 = 2.5 GeV 2 ) are plotted on Fig.3. Two sum rules give the same value of A^ at 
a = 0.55 GeV 3 . The 10% variation of W 2 does not change this result. The M 2 -dependence 
of X 2 N , determined from (A.l) and (A. 2) at these values of fitting parameters is shown on 
Fig.4. As is seen, Afy found from two sum rules agree with one another with accuracy 
~ 3% andjheir deviation from constant is less than 5%. The mean value of \ 2 N can be 
chosen as X% — 2.1 GeV 6 {n = 1 GeV). 
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Figure Captions 



Fig. 1. The independence of £ + C M 2 at / 2 = 3 x 1(T 2 GeV 2 , eq.25, 
g\ + C 8 M 2 , and g A - 1 + C^M 2 , eq.29. 

Fig. 2. £ (solid line, left ordinate axis) and \^x z , eq.(27), (crossed line, 
right ordinate axis), as a functions of /g. 

Fig. 3. The values of as functions of a determined from the sum rules 
(A.l) - solid line and (A. 2) - crossed line. 

Fig. 4. The M 2 - dependence of Afy found from the sum rules (A.l) - solid 
line and (A. 2) - crossed line. 
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